Joint superexchange— Jahn- Teller mechanism for A-type antiferromagnetism in LaMnC>3 
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We propose a mechanism for A-type antiferromagnetism in orthorombic LaMn03, compatible 
with the large Jahn- Teller splitting inferred from structural data. Orbital ordering resulting from 
Jahn- Teller distortions effectively leads to A-type ordering (antiferromagnetic in the c axis and 
ferromagnetic in the ab plane) provided the in-plane distorsion Q2 is large enough, a condition 
generally fulfilled in existing data. 
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Stoichiometric LaMn03 (LMO) is known |jj to be 
an A-type antiferromagnetic insulator (A-AFMI), where 
ferromagnetically ordered Mn02 planes (in the xy di- 
rection) have staggered magnetization along the z axis. 
Upon increasing the temperature a paramagnetic insulat- 
ing phase (PMI) is reached. On the other hand sufficient 
hole doping (e.g. by substituting La with Sr or Ca) gives 
rise via the so-called double-exchange hopping mecha- 
nism |^|,[| to a low-temperature ferromagnetic metallic 
phase (FMM) turning into a PMI phase at higher temper- 
ature. Not only magnetism determines the main physical 
properties, in fact both theoretical and experimental 
H evidences emphasize the relevance of electron-lattice 
coupling. Charge and orbital ordering also occur, fur- 
ther showing the competition between various physical 
mechanisms. Notice that the crucial role of spin and lat- 
tice coupling was repeatedly emphasized to account for 
the properties of the FMM phase and the FMM-PMI 
transition at finite doping, as well as charge-ordering 
phenomena. However, this liaison regarded the double- 
exchange mechanism for charge tranport, being dynam- 
ically dressed by lattice degrees of freedom No 
emphasis was put on the role of static cooperative Jahn- 
Teller (JT) deformations in stabilizing specific magnetic 
structures in the AFMI phase. 

Here we investigate an approach to the stoichiometric 
phase of LMO showing that the layered antiferromag- 
netic structure may result from the interplay between 
superexchange and JT couplings. Our analysis is alter- 
native to the more qualitative one based on the semico- 
valent exchange mechanism M and is complementary to 
the superexchange mechanism investigated by Kugel and 
Khomskii (KK) || for perovskites with JT ions. This lat- 
ter analysis (see also |)| ) focused on the interplay between 
magnetic and orbital ordering within the two e g orbital 
manifold, assumed degenerate. In particular, basic ingre- 
dients were a strong local electron-electron repulsion U, 
the Hund coupling Jh between electrons on the two e g 
orbitals, and the orbital mixing (described by a mixing 
angle 8) due to JT distorsion. In the approach of Ref . || , 
only e g (spin and orbital) degrees of freedom were consid- 
ered, and the spin and orbital order were self-consistently 



determined to lowest order in Jh/U. The e g level degen- 
eracy was lifted by superexchange but the JT splitting 
resulting from the lattice distorsions induced by orbital 
ordering was not explicitly considered, Only a correction 
due to a small local JT anharmonicity was introduced. 
This point of view, considering magnetic exchange to be 
the main cause of orbital mixing/ordering, but neglect- 
ing the orbital splitting resulting from JT effect, might 
be questioned in LMO where strong JT distorsions arise. 
Moreover, magnetic exchange interactions are somewhat 
modest, e.g. Jaf ~ 0.58 meV and Jf ~ 0.83 meV > Jaf 
from inelastic neutron scattering experiments [ [Tof , while 
the KK theory results in Jp « (Jh /U)Jaf, that is Jp 
much smaller than Jaf- Moreover, as pointed out by 
KK, the observed orthorhombic distorsion with c < a 
is not expected from considering solely orbital ordering. 
Therefore, it is not obvious, considering the actual distor- 
sion, that the magnetic A-phase is still the most stable. 

According to the experimental evidences for the rele- 
vance of the t2 g (spin) degrees of freedom, e.g. in the 
double exchange hopping processes, and for a strong JT 
coupling, we propose to reconsider the problem. We take 
properly into account the Hund coupling between e g and 
ti g electrons, and assume, contrarily to KK, that the JT 
splitting is much larger than the exchange energy. This 
makes the JT effect the driving mechanism for orbital 
ordering, which in turn controls the magnetic interac- 
tions. Hereafter we assume some staggered orbital order 
(i.e. we fix 8 on the two sublattices) as determined by 
strong JT distorsions, without attempting to calculate 
them, since they can be extracted from crystallographic 
data. Given the orbital order and the related (relatively 
large) e g orbital splitting e, we calculate at any order 
in e and Jh/U the superexchange interactions. Finally, 
we determine the parameter ranges which are compatible 
with the observed A-type AF phase. 

For the sake of simplicity we disregard the oxygen sites 
in the perovskite structure, thus focusing on a single- 
site model. On each site two Manganese orbitals, the 
d x 2_ y 2 (x) and the d^ z 2_ r 2 (z), are available. The real 
lattice structure is effectively taken into account via the 
sign and the magnitude of the intersite hopping along the 
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x — y — z directions. Specifically, we notice that, for a 
standard choice of phases, the hopping between the d x 
and the d z orbitals on Mn are given by 

— V3i along x 
Vst, along y 
. 
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Throughout this paper the apex indicates hoppings in the 
z direction. The (static) JT distortions |ll]] mix the x and 
z orbitals [E| into a (lower) and b (upper) orbitals with e g 
symmetry split by an energy e = 2g\/Q\ + Q\ where g 
is an electron-lattice coupling constant. Specifically, us- 
ing the same notations of Ref. j^] the uniform distortion 
Qi (corresponding to a uniform variation of the lattice 
parameter along the z direction) couples to the x — z den- 
sity difference n x — n z . In systems like LMO, where the 
lattice spacing in the z direction is shorter that the (aver- 
age) spacing in the xy plane, the d z orbitals overlap more 
and are pushed at higher energy by the Coulomb inter- 
action. Then Qi is negative. As pointed out by KK, this 
disfavours hopping in the z-direction, thus making the 
JT effect compete with the A-type superexchange which 
alone would imply the opposite distortion. The reverse is 
true in systems like KCUF3, where the lattice parameter 
along z is larger than in the xy plane and Qi is positive. 
On the other hand, the distortion Q2 (opposite on the 
two sublattices of the xy-planes) corresponds to an alter- 
nate contraction and dilation of the Mn-0 bonds on the 
xy plane and mixes the x and z components of the a, b 
orbitals 



I a) = cos(0/2)|a;) ± sin(0/2)|z) 
1 6) = sin(0/2)|:c) =f cos(0/2)|z) 



(2) 



where tan(#) = Q2/Q1 and the upper (lower) sign is for 
sites on sublattice 1 (2) of the xy-planes. Accordingly 
the hoppings between the a and b orbitals of neighbour- 
ing sites can straightforwardly be obtained via Eqs. (|l|). 

t aa = -< (1 + 2 cos 61); (3) 
t ab = -t(±V3 + 2sinl9); t ba = -t (+V3 + 2 sine) ; (4) 
C = 2t(l-cos0); t' ab = t' ba = -2tsm6 (5) 

the upper (lower) sign is for planar hopping in the x (y) 
direction. 

As customarily done, we assume that the Hund cou- 
pling between a — 1/2 e g electrons and the 5 = 3/2 spin 
of the t2 g electrons is so large that the initial and final 
states always have maximal total spin St — 2. Moreover, 
and most importantly, we also consider a large local re- 
pulsion (~ U) between electrons on Mn sites forbidding 
to two electrons to reside on the two e g levels of the same 
site. Then we work on a reduced Hilbert space with only 
Mn 3+ initial and final states [fill 



We then carry out a perturbative calculation of both 
FM and AFM magnetic couplings between sites 1 and 2 
by considering second order hopping processes from and 
to the ground state configuration with one electron per 
Mn site occupying the lower a orbital. We thus neglect 
the exchange-induced mixing of e g orbitals, considered 
by KK. Notice that this last assumption relies on the JT 
splitting e being substantially larger than both the tem- 
perature and the superexchange scale ~ t 2 /U. Due to 
the condition St = 2, each site i = 1, 2 is five times de- 
generate, |2, m)i with m = —2,-1,. ..,2. The first step 
consists in forming two-site states with given total spin 
J = 0, 4 from the 25 basis states |2, m)i® |2, 771)2- The 
suitable Clebsch-Gordan coefficient are easily obtained. 
Within each J subspace (the hopping processes conserve 
the total spin), the hopping perturbation 



H t = - 



^ ^ taa' 
(T;a,a'— a,b 



is applied twice to obtain the (J, M\Hf\J, M) matrix el- 
ements. The double-hopping processes are of two types: 
aa and ab depending on whether the initial a electron 
hops on the neighboring a or b orbitals. Accordingly 
there are two superexchange channels, leading to cou- 
plings constants J aa and J a t,. The easiest to be calcu- 
lated is J aa since the Pauli principle forces the two initial 
a electrons to have opposite spins. As a consequence only 
one intermediate virtual state per J channel is allowed, 
with one empty and one doubly occupied a orbital. Both 
the empty and the doubly occupied orbitals cost an en- 
ergy 3J^/4, where J' H is the Hund coupling between e g 
and t2 g orbitals (the Hund energy is set to zero in the 
Mn 3+ ground state configuration). The doubly occupied 
a orbital has an additional energy cost U. All the inter- 
mediate states have then an energy Ey = U + (3/2) J' H 
above the ground state energy E$ = 0. As it is standard, 
the perturbative energy gain can be compared with the 
energies of | J, M) states as given by the effective Hcisen- 
berg model for the aa channel 

H aa — Jaa (Si ■ S2 + C) 

with C being a constant energy shift to be determined. 
The direct comparison provides, besides C = — 4, 



J aa 



1 



i, 



AU + ('S/2)J' H 



(G) 



This is an effective AFM coupling between electrons on 
the xy-plane. A similar expression is obtained for the in- 
terplane coupling (i.e., in the z direction) J' aa , provided 
t aa is replaced by t' aa in Eq. (|). 

The calculation for J a b is slightly more complicated, 
since the hopping electron can now give rise on the dou- 
bly occupied site to both a, St — 5/2 or a St = 3/2 
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state, thus increasing the number of virtual states. How- 
ever, the same procedure illustrated above yields 



j ab 
Jab = 77" 



' ab 



40 \U' + e-J H /2 



U' + e- J H /2 + W H /2 U> + e + J H /2 + 3J' H /2 

(7) 

where U' is the local Coulomb repulsion between elec- 
trons on the a and b orbitals and 



f z = 



_(* 2 «6 



tla) 



= t 2 (3 + 4 sin 2 (0)) . 



(8) 



Again the analogous coupling in the z direction can be 
obtained by replacing t ab and tt, a with the corresponding 
primed quantities of Eq. (|^) . It can easily be seen that 
this coupling is ferromagnetic and vanishes when Jh and 
J' H are both zero. Notice also that, since J a b arises from 
virtual hopping i — > j — > i and j — ► i — ► j and since the 
a and b orbital combinations are reversed on neighboring 
sites, the t 2 b combination appears, which is the same in 
the x and y directions [cf. instead (|])]. Thus for each 
crystalline direction one can write the effective Heisen- 



berg couplings J = J aa + Jab and J' = J' aa + J' at 



the 



xy planes and z direction respectively. Then the question 
arises concerning the parameter ranges such that the ob- 
served A-type AF is realized. In this case the coupling 
must be dominantly ferromagnetic in the xy planes and 
dominantly AF in the z direction. To this purpose we 
rewrite the J's in the following more compact way 



Jaa — t a a/D 



1 1 in Jab — t ab / Dab 



(9) 



The condition that the xy planes are ferromagnetically 
coupled is written as \J a b\ > Jaa, i-e. 



a 2 = 



t 2 

b ab 

t 2 
'a a 



> D. 



(10) 



where D — At the same time, the condition for AF 
coupling in the z direction is expressed by \J' ab \ < J' aa , 
that is 



t'. 



< D. 



(11) 



Now, both a and a' (i.e. t aa , t a b, t' aa , and t' ab ) are frui- 
tions of 9, or of the JT ratio Q2/Qi- Plotting a 2 and a' 2 
as a function of Q2/IQ1I f° r Qi < (the relevant case for 
LMO), one obtains the curves in Fig. 1. Since from the 
inequalities dlTI) and (fll), one can deduce the condition 



a 2 > D > a' 2 , 



(12) 



the value of D should be below the solid curve and above 
the dashed ones. The assumed orbital order becomes 



compatible with A-AFMI order for Q2/IQ1I 2.6, where 
the relation (|TJ|) can be satisfied. Notice instead that for 
Qi > the curve a' 2 is always smaller than a 2 , so that 
no restriction on the Q2/Q1 ratio is needed to fulfill the 
condition (|l2]). This is in agreement with the coopera- 
tion between superexchange and JT effect, found in this 
case, see for example the case of KCuF 3 H. 
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Fig. 1.: a 2 (solid line) and a' 2 (dashed line) as a function of 
the JT ratio Q2/IQ1I for negative Qi. The dotted line rep- 
resents a' 2 with t' xz = 0.05t. The vertical lines indicate the 
experimental values Q2/IQ1I = 2.84, 6.69 (see text) 



More precisely we find that 



J 
J 7 



t 2 D 



a 



t'„ 



D - a'- 



(13) 



Notice that a substantial amount of Q 2 /\Qi\ JT dis- 
tortion is needed to leave the possility open for the con- 
dition ( |T2| ) to be fulfilled. Using standard results 111 ] 
connecting the structural parameters with the Q2/IQ1I 
ratio (Q2/IQ3I m the notation of Ref. p2[ ) we estimated 
Q2/IQ1I ~ 1-97 and 6.69 for the parameters of the or- 
thorhombic structures reported in Table V of Ref. p[ , 
while a value 2.84 is obtained from the data related to 
the crystal where the abovementioned values of Jaf and 
Jf have been measured pfj . Although not all values are 
compatible with the condition Q2/IQ1I ~ 2.6, the large 
variability of the resulting estimates for Q 2 /\Qi \ indicates 
that this parameter may vary from a crystal to another 
so that the condition Q2/IQ1I ?i 2.6 is quite reasonable. 
Especially, we now show that the last one yields values 
consistent with the observed J's. 

Let us determine the (hopefully realistic) values of U , 
U', Jh, J'h an d e providing a D ratio in the needed 
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range (|l2|). To this purpose we consider typical val- 
ues of U = 10 eV, U' = U - 2J H = 6 - 10 eV, and 
e = 0.1 — 0.5 eV, and we obtain the ratio D as a function 
of Jh/U. For simplicity the ratio Jh/J'h is taken equal 
to 1. Then, from Fig. 1, we consider typical ranges for 
D between, e.g., a' 2 rs 1.4 and a 2 rs 4., roughly corre- 
sponding to Q2/IQ1I ~ 6.69, and a' 2 « 2. and a 2 w 2.4, 
roughly corresponding to Qs/|Qi| ~ 2.84. Then we look 
for what range of Jh/U the conditions 4 > D > 1.4 or 
2.4 > D > 2. are realized. In the first case the quite 
reasonable range 0.1 < Jh/U < 0.35 is found, while in 
the second case 0.1 < Jh/U < 0.15. The observed values 
J = 0.83 meV (ferro) and J' = 0.58 meV (antiferro) 
are obtained with Jh — J'h — 1-25 eV and t = 0.34eV in 
the case of Q2/IQ1I = 2.84. This value of t is fairly large, 
but it can be substantially decreased by taking larger 
values of Q2/\Qi\- Notice that, contrarily to the cal- 
culations of KK, the planar ferromagnetic coupling can 
easily be larger than the z-axis antiferromagnetic one. 
This comes from the ratio t 2 a /t' aa 2 in (|p]) where, due 
to the orbital order stabilized by the JT distortion, hop- 
ping in the plane directions is enhanced with respect to 
the z- axis. 

In conclusion, we have illustrated the possibility of 
an alternative mechanism for the layered antiferromag- 
netism of the stoichimetric LMO compound. In partic- 
ular, we showed that the superexchange mechanism, to- 
gether with strong JT planar distortions, can be respon- 
sible for the specific A-type magnetic structure. Like 
in the analysis of Ref. 13], the sinergetic effect of both 
magnetic superexchange and orbital ordering is a crucial 
ingredient. However, contrary to the assumption of Ref. 

in the present scenario, we assumed a given orbital 
ordering strongly lifting the degeneracy of the e g orbitals 
(e S> J a b, J aa )- A relevant role here is also played by the 
tgg spin degrees of freedom, as seen from the expressions 
(||) and (Q), which depend rather strongly on J' H . 

Despite the basic differences between our scheme and 
the proposal of Ref. [|| , we find some similarities in the 
overall result. Specifically we find that for systems with 
positive Qx, like, e.g., KCUF3, no restriction is needed on 
the Q2/Q1 ratio to make the orbital ordering compatible 
with A-AFMI magnetic structure. This is not the case 
for negative Qi, where the above discussed conditions 
have to be imposed on Q1/IQ2I and, consequently on the 
values of Jh/U. Similar values of the Jh/U were found 
in |. 

Within the present model, to be more realistic one 
should also take into account the existence of tilt distor- 
tions and consider the effects of non vanishing transfer 
integrals t' xx and t' xz . In fact these hopping constants are 
strictly zero only for lattices without tilting of the MnOg 
octaedra around axes on the xy plane. We have found 
that positive t' xz and negative t' xx favour A-type ordering. 
Another contribution to the problem is the antiferromag- 
netic exchange originating from t2 g electrons. Finally, 



a full calculation involving all superexchange processes 
should be feasible. 

We stress again that the mechanism proposed here is 
based on the tight interplay of lattice and electronic de- 
grees of freedom existing even in the undoped LMO sys- 
tem, and that it is aimed to correlate the magnetic or- 
dering with the JT distortions. Justifying the values of 
Qi and Q2 from microscopic grounds is beyond the scope 
of the present paper. Some mixing between lattice and 
spin dynamics and isotopic or pressure dependence of the 
spin- wave velocity are expected to be rather natural con- 
sequences of the proposed scenario. 



ACKNOWLEDGMENTS 

M.G. acknowledges financial support of the Istituto 
Nazionale di Fisica della Materia, Progetto di Ricerca 
Avanzata 1996, as well as from Universite Joseph Fourier 
and LEPES, CNRS in Grenoble where part of this work 
was carried out. D.F. acknowledges hospitality and sup- 
port from the Dipartimento di Fisica of Universita di 
Roma "La Sapienza" . 



[1] E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 

[2] C. Zener, Phys. Rev. 82, 403 (1951). 
[3] P. W. Anderson and H. Hasegawa, Phys. Rev. 100, 675 
(1955). 

[4] A. J. Millis, P. B. Littlewood, and B. I. Shraiman, Phys. 
Rev. Lett. 74, 5144 (1995); A. J. Millis, R. Mueller, and 
B. I. Shraiman, Phys. Rev. B 54, 5405 (1996) and refer- 
ences therein. 

[5] J. Zang, A. R. Bishop, and H. Roder, Phys. rev. B 53, 
R8840 (1996). 

[6] For a recent review see, e.g., A. R. Bishop and H. Roder, 

|cond-mat 970314S 
[7] J. B. Goodenough, Phys. Rev. B 100, 564 (1955). 
[8] K. I. Kugel and D. I. Khomskii, Sov. Phys.-JETP, 37, 

725 (1973). 

[9] M. Cyrot, C. Lyon-Caen, J. Physique 36, 253 (1975). 

[10] F. Moussa, et al., Phys. Rev. B 54, 15149 (1996). 

[11] Y. Yamada et al, Phys. Rev. Lett. 77, 904 (1996). 

[12] J. B. Goodenough, Magnetism and the Chemical Bond, 
(Ch. Ill, Par. IF), Ed. Interscience Publishers 1963, John 
Wiley and Sons, New York - London. 

[13] This is consistent with the observation that, in the ab- 
sence of strong-correlation effects, the ferromagnetically 
ordered Mn02 planes at quarter filling (one electron and 
two e g orbitals per site) should be metallic (unless an 
exceedingly large JT splitting e of the order of the band- 
width or more is assumed). 

[14] J. F. Mitchell, et at, Phys. Rev. B 54, 6172 (1996). 



4 



